arXiv:hep-ph/0112139v2 12 Apr 2002 


Production of intermediate-mass dileptons in relativistic heavy ion collisions 


loulia Kvasnikova and Charles Gale 
Department of Physics, McGill University, 

3600 University Street, Montreal, QC, H3A 2T8 Canada 

Dinesh Kumar Srivastava* 

Department of Physics, Duke University, Durham, NC 27705 
(February 1, 2008) 

The production of intermediate mass dileptons in ultrarelativistic nuclear collisions at SPS en¬ 
ergies is studied. The acceptance and detector resolution inherent to measurements by the NASO 
experimental collaboration are accurately modeled. The measured centrality dependence of the 
intermediate mass lepton pair excess is also addressed. 


I. INTRODUCTION 

The study of relativistic heavy ion collisions derives its unique importance from the opportunity of producing hot 
and dense strongly interacting matter. More specifically, one of the hopes is the creation and study of a quark- 
gluon plasma (QGP), a prediction of QCD Q, which filled the nascent Universe roughly a microsecond after the Big 
Bang. This field of endeavour currently represents a cutting edge of research in subatomic physics and has generated 
tremendous activity, both in theory and in experiment . Good penetrating probes of in-medium phenomena in heavy 
ion collisions are electromagnetic observables, owing to their small rescattering cross sections and to their privileged 
coupling to vector mesons |^. For collisions at ultrarelativistic energies, the lepton pair invariant mass (M) spectrum 
can roughly be divided in three regions: low mass (M < m^), high mass (M > and intermediate mass. The 

high mass region is Drell-Yan-dominated at the CERN SPS, and this contribution is in principle calculable within 
the realm of perturbative QCD. In connection with the low mass sector in heavy ion collisions at similar energies, the 
measurements of the CERES collaboration Q have been seminal in building the case for an observation of in-medium 
effects on the vector meson spectral densities |^. There, the experimental collaboration has identified an excess over 
the sources that were known to be sufficient for the understanding of the spectra in pA collisions. This modification 
of the vector meson properties in matter represents an exciting discovery and is a genuine consequence of many-body 
physics. 

In the intermediate invariant mass region, an excess of dimuons over the sources expected from pA measurements 
has also been identified experimentally by the Helios/3 and NA50 Q| collaborations. Some excitement and interest 
have also been associated with this observation as this mass region is the one originally predicted to contain a plasma 
signal Q. Theoretical models have been devised to understand this dimuon excess. Those include rescattering of 
the open charm mesons with the surrounding hadronic constituents [^. This approach can not however reproduce 
the dimuon mass and transverse momentum spectra for central nuclear collisions 0. An overall charm production 
enhancement has also been invoked as a prospective interpretation [^. Whereas this scenario raises an interesting 
possibility, it has to be reconciled quantitatively with the observation of J/^/> suppression. This issue can however 
be settled by a direct measurement |ll[. An earlier theoretical study of the Helios/3 measurements has attributed 
the excess to secondary hadronic reaction processes |l^. This result has some consequences that are important to 
put on a firmer basis, by comparing with other newer data. In this strain, more recent NA50 results have also been 
interpreted by thermal models, with some quark-gluon content 00 ■ 

Our goal in this paper is to consider the NA50 intermediate mass dimuon results in the light of a hydrodynamic 
modeling of the nuclear dynamics, of a detailed analysis of the hadronic dilepton rates, and of a precise simulation 
of the detector acceptance and resolution. The issue of centrality dependence of the signal is also addressed. The 
next section briefly describes the hydrodynamic approach utilized, along with its reproduction of hadron spectra. The 
dilepton rate calculation is then described and the numerical modeling of the detector cuts and resolution is introduced. 
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After the dilepton spectra in invariant mass and transverse momentum have been discussed, the centrality dependence 
is considered. The paper closes with a summary and a conclusion. 


II. HYDRODYNAMIC EVOLUTION AND HADRON SPECTRA 


It is well known that the final yield of the thermal dileptons has to be obtained by an integration over the space-time 
history of evolution of the interacting system. We shall assume that a thermally and chemically equilibrated quark- 
gluon plasma is produced in such collisions at a time tq, which may be estimated from the condition of isentropic 
expansion of the plasma 11]: 


27r4 I dN 

4 ^ 


( 2 . 1 ) 


where dN/dy is the particle rapidity density for the collision, and a = 42.257r^/90 for a plasma of massless u, d, s 
quarks, and gluons. The number of flavours is taken to be Nf Ki 2.5 to account for the mass of strange quarks. We 
see that once the transverse area of the colliding system Ax is known, along with dN/dy, the above relation uniquely 
relates Tq to tq. While one may get arbitrary large value of Tq by choosing a smaller formation time of the plasma, 
the uncertainty principle provides a lower physical limit on tq ~ l/STg |l^ . 

While considering Pb -I- Pb collisions at SPS energies, we have taken the average particle rapidity density as 750 
for the 10% most central Pb -I- Pb collisions at the CERN SPS energy, as measured. We estimate the average number 
of participants for the corresponding range of impact parameters (0 < b < 4.5 fm) as about 380, compared to the 
maximum of 416 for a head-on collision. We thus use a mass number of 190 to get the radius of the transverse area of 
the colliding system and neglect its deviations from azimuthal symmetry, for simplicity. As the deviation measured 
in terms of the number of participants is marginal (< 9%), we expect the error involved also to be small. We also 
recall that the azimuthal flow is minimal for near central collisions. 

The plasma is then assumed to undergo a boost-invariant longitudinal expansion and an azimuthally symmetric 
radial expansion, with a transition to a hot hadronic gas consisting of all hadrons having M< 2.5 GeV, in a thermal an 
chemical equilibrium at a transition temperature T^. This makes for a rich equation of state. Once all parton matter 
is converted into hadronic matter, we assume the hot hadronic system to continue to expand until it undergoes a 
freeze-out at some temperature Tp. During this evolution, the speed of sound of the matter is consistently calculated 
at every temperature to be used in the equation of state needed for solving the hydrodynamic equations JT^ . 

It is important to realize that the (pressure) gradients in the system cause and control the radial expansion of 
the plasma, which in turn leads to a more rapid cooling and to the development of a transverse velocity for the 
fluid. Support for these ideas comes from intensity interferometry and from measured particle spectra. It is thereby 
necessary to have a proper initial energy density profile as it affects the hydrodynamic developments by introducing 
additional gradients. We assume it to follow the so-called “wounded-nucleon” distribution, which for central collision 
of identical nuclei leads to: 


e(To,r) oc / p{\/r^~+^)dz (2.2) 

J —OO 

where p is the (Woods-Saxon) distribution of nucleons and r is the transverse distance. This is prompted by the 
experimental observation that transverse energy deposited in these collisions scales with the number of participants. 
The normalization is then determined through numerical integration such that 


41t Co = J 27rre(r) dr. (2.3) 

A proper profile is also needed to get a quantitative description of the dilepton yields and through this to infer the 
initial conditions. It is clear that the energy densities at larger r would be smaller, implying lower temperatures and 
a shorter life-time till the fluid there undergoes freeze-out. Now recall that the four-volume d'^x = t dr r dr dp d(j) and 
thus the contribution of a given fluid element rises linearly with its transverse distance, which further “magnifies” 
the contribution from large r. One can immediately see that a uniform energy density profile often employed in the 
literature would lead to an erroneously large contribution from large r. This fact plays an important role in arriving 
at the required initial conditions. 

We now have all the ingredients to solve the hydrodynamic equations. This is done through the procedures described 

in Refs. 1T1,|T|. 
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A. Initial conditions at the SPS 


On the basis of the rapidity density etc. observed at SPS energies in central collisions involving lead nuclei, it was 
shown recently that a good empirical description of the single photon data measured by the WA98 collaboration 
was obtained when the plasma formation time was taken to be tq = 0.2 fm/c. This corresponded to an average initial 
temperature of about 330 MeV. In the present work, we shall use these values but also explore the consequences of 
varying tq and Tq, keeping the resulting dN/dy fixed. 

The phase transition is assumed to take place at T = 180 MeV and the freeze-out at 120 MeV. This value of 
the critical temperature is motivated by lattice QCD results which give values of about 170 - 190 MeV [|], and the 
thermal model analyses of hadronic ratios which suggest that the chemical freeze-out in such collisions takes place at 
about 170 MeV (recent analyses yield a value of 158.1 ± 3.2 MeV |2^] for the chemical freeze-out temperature). The 
phase transition should necessarily take place at a higher temperature. 


B. Hadronic Spectra 


As a first step we look at the spectra of several hadrons for the central collision of lead nuclei measured at the SPS. 



FIG. 1. Hadronic spectra in central collision of lead nnclei at SPS energy. The formation time is assnmed to be 0.2 fm/c. 
The data are for pions, protons, kaons, and lambdas -|- sigmas, from top to bottom, respectively. 



FIG. 2. Hadronic spectra in central collision of lead nuclei at SPS energy, 
fm/c respectively. 


The formation time is taken as 0.2 fm/c and 1 
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In Fig. 1^ we show our results when the formation time is take as tq = 0.2 fm/c as in Ref. . The data are from Refs. 

Even though it is known for quite some time that a reasonable variation of the initial temperature and the 
formation time, keeping the corresponding dN/dy fixed, affects the flow only marginally, we study the consequences of 
increasing tq to 1 fm/c in Fig. |[ We see that the hadronic spectra shapes are relatively insensitive to such variations 
in the initial conditions, as the flow takes some time to develop. In all cases, the agreement is quite good. This fact 
lends some credibility to the dynamics advocated here. 


III. DILEPTON SPECTRA AT SPS ENERGY 


A. Dilepton emission rates 


For the dilepton emission from the deconfined sector of QCD, we have assumed that the dominating process is 
qq —> The emission rates that correspond to this process have been calculated many times. An example of 

their quantitative contribution can be found in Ref. [ p^ . Even though those rates are the Born terms in the parton 
sector, a very recent lattice calculation has found that those were surprisingly accurate |25|. 

In the evaluation of hadronic rates most calculations have relied on effective Lagrangian techniques, where the 
available parameters are fitted to empirical measurements. An example of the application of this line of thought to 
the low mass dilepton sector can be found in Refs. |^,^. A problem which arises when one tries to extend those 
techniques to the intermediate mass region is the appearance of off-shell effects. Indeed, different approaches that 
agree in the soft sector can yield widely different results in higher invariant mass extrapolations P7[ . Fortunately, 
constraints on the hadronic processes can be obtained through the wealth of data of the type e+e” ^ hadrons [ p^ . 
In addition, those measurements cover exactly the same invariant mass range as the one which concerns this work. 
Those have been used, together with r-decay data, to construct the vector and axial-vector spectral densities [p^ . 
Similarly, the intermediate invariant mass initial state e+e“ data has been analyzed specifically in the most important 
exclusive channels. This information can then be used to construct rates for hadrons ^ e+e“ [|^. It is this procedure 
that is followed in the current work. The contributing channels producing lepton pairs in the invariant mass range I 
GeV < M < 3 GeV have been found to correspond to the initial states: tttt, irp, ttw, yp, pp, ttoi, KK, KK* + c.c. 
[p^ , p0| . A detailed channel-by-channel discussion is too long to be had here, but a quantitative assessment of the net 
dilepton rate from a gas of mesons at T = 150 MeV is shown in Fig. ^ Also shown is the rate extracted from the 
spectral function evaluation In principle, the latter contains all the hadronic sources. The correspondence 

between those two results is a confirmation that the important channels have indeed been identified in the hadronic 
scenario. 


> 

0 ) 





FIG. 3. Net dilepton production rate from a gas of mesons at a temperature of T = 150 MeV, as a function of dilepton 
invariant mass. The full curve is the sum of the hadronic channels discussed in the text and in the references. The data points 
follow from an extraction of spectral densities from e'''e~ scattering experiments @,13 • 
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Using the rates of dilepton production from quark and hadronic matter, it is straightforward to calculate the spectra 
for dileptons, say for the central collision of lead nuclei at SPS energy. We add the contributions of the quark matter 
from the QGP phase and the quark matter part of the mixed phase and call it QM. Similarly we add the contributions 
of the hadronic matter part of the mixed phase and the hadronic phase and call it HM. The contribution of the Drell- 
Yan process is obtained by a scaling of nucleon-nucleon estimates to the case of nncleus-nucleus collisions. We get the 
results shown in Fig. no acceptance or resolution corrections are implemented. It is clear that in principle, if the 
formation time is small (ie., the initial temperature is large) then there could be a substantial contribution of the 
dileptons having their origin in the quark matter, in the intermediate mass window. Note that the HM contribution 
is only very marginally affected by variations in the tq. We shall proceed to model the effect of detector resolution 
and to compare with experimental measurements. 



M (GeV) 


FIG. 4. The net yield of lepton pairs from the quark gluon plasma, Drell-Yan and the hadronic sources are shown separately, 
before detector and acceptance corrections. The QGP signal for different formation times tq are shown (solid lines), earlier 
time on top, later times on the bottom. The Drell-Yan contribution is the dashed line, and the hadronic matter contribution 
the dashed-dotted curve. The yield from correlated open charm decay is not shown here. That source’s contribution to dN/dM 
is a decaying exponential in the mass range under scrutiny and its dependence on various parameter assumptions is discussed 
in detail in [M. 


B. Detector acceptance and dilepton spectra 

It is vital to account for the finite acceptance of the detectors and for their resolution when comparing the results 
of theoretical calculations with measured experimental data. Those effects are important in the NASO experiment 
0-. One approach to this problem is to model approximately and analytically the acceptance. While this can be 
readily implemented 0,^, a legitimate doubt can subsist about the accuracy of the experimental representation, 
especially in regions where edge effects might be important. In order to circumvent this problem we use here for the 
first time a numerical subroutine developed to reproduce the NASO acceptance cuts and finite resolution effects in the 
measurements of dimuon pairs in Pb -|- Pb collisions at the CERN SPS [0| . 

We compute the invariant mass distribution of the dileptons in our hydrodynamic approach and then we run our 
pairs through the numerical detector simulation. The normalization is determined by a fit to the Drell-Yan data using 
the MRSA parton distribution functions as in the NASO analysis. For getting the px distribution we supplement our 
dN/dM"^ estimates for the Drell-Yan with a gaussian distribution in px as in Ref. [0 which very closely reproduced 
the estimates obtained by the NASO collaboration. The result is shown in Fig. |[ We also compute the px distribution 
for the muon pairs, in the invariant mass window l.S < M < 2.S GeV. This spectrum is broken down in the same 
sources as previously and is shown in Fig. 
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FIG. 5. We show the calculated dimuon invariant mass spectrum after correcting for detector acceptance and resolution. 
The data is from the NASO collaboration |Q. The Drell-Yan and thermal contributions are shown separately, as well as those 
coming from correlated charm decay and direct decays of the J/'ip and ip' . 

At this point it is appropriate to consider the following question: which initial temperature is demanded by the 
intermediate mass dilepton data? A critical and quantitative assessment of this issue can be obtained by examining 
a linear plot of the lepton pair spectrum in the mass region under scrutiny. This is show in Fig. From this figure, 
it appears that the best fit is provided by tq = 0.2 fm/c, and that the second best (less than two standard deviations 
away for most of the data points) belongs to tq — 0.4 fm/c. In terms of initial temperature those two values of the 
formation time correspond to Tq ss 330 and 265 MeV respectively. A conservative and reasonable point of view is that 
it is probably not fair in such a challenging and complex environment as that of ultrarelativistic heavy ion collisions 
to ask for an agreement that is better than two standard deviations, considering the inherent uncertainties. This 
circumvents the issue of initial temperature determination. Also, there are questions that have to do with the specific 
dynamics used here to model the nuclear collisions. However in this respect, the hydrodynamic model used in this 
work is solved self-consistently and its parameters are constrained by the hadronic data as well. In connection with 
the important discussion on plasma signals the quark matter contribution is ~ 23% for tq = 0.2 fm/c, and « 19% for 
To = 0.4 fm/c, around the lepton pair invariant mass of 1.5 GeV. 



FIG. 6. We show the dimuon transverse momentum spectrum after accounting for detector effects. The origin of the data 
and of the different sources are as in Fig. 
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FIG. 7. A linear plot of the net dilepton spectrum in the intermediate mass region. The three solid curves correspond to 
formation times tq = 0.2, 0.4, 0.6 fm/c, from top to bottom, respectively. The data are from Q. The thermal contributions 
and the contribution from hard processes are shown separately. 


C. Centrality dependence 

Our discussion so far was concerned only with the high multiplicity bin or, with mostly central collisions. To extend 
the hydrodynamic model to non-central events and to properly treat the azimuthal anisotropy is not a simple task. 
However, one can get a fair estimate of the centrality dependence by ignoring the broken azimuthal symmetry and 
by approximating the region of nuclear overlap by a circle of radius R « 1.2(A^part/2)^'^^, where iVpart is the number 
of participants ||34j| . Assuming the uncertainty relationship, tq = l/3To, one can track the changing multiplicities 
without introducing additional parameters. The results of this exercise are shown in Fig. for S -I- U and Pb -I- 
Pb collisions. Also shown are the contributions from the “enhanced production of charm” estimated by the NA50 
collaboration, normalized to the dileptons estimated by us. We estimate the contribution of the “excess charm” ANch 
from 


ANch oc {E - l)Ncou 


(3.1) 


where Ncoii is the av erag e number of collisions for the given centrality and E is the ‘enhancement’ factor given by the 
NA50 measurement [Q. The constant of proportionality is identical for all the bins for a given system. It is seen 
that the agreement with the measured data is quite good, and our approach gives a fair description of the centrality 
dependence of the excess dilepton measurement. 



FIG. 8. Gentrality dependence. The data represent the “extra charm” yield (as categorized by the NASO collaboration Q) 
needed to describe the intermediate mass region dimuon data. The full curve is from our thermal sources. 
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IV. SUMMARY AND CONCLUSION 


In the light of the findings reported on in this work, it does appear that the case of attributing the intermediate 
mass dilepton excess to an enhanced charm production rests on evidence that is no longer very compelling. The lepton 
pairs coming from a thermal hadronic medium comprise an important part of the signal, in agreement with earlier 
analyses of the Helios/3 measurements |^2|, and with those of more recent studies |^,Q. It is satisfying to notice 
that the behaviour of the centrality dependence is also reproduced by the modeling and the dynamics used here. We 
have also considered the constraints on the initial temperature set by the data in the framework of our approach. 
There is a QGP component in the treatment used in this work. Unfortunately, its numerical value is not important 
enough to provide irrefutable evidence. The future of this field is bright: RHIC electromagnetic measurements should 
soon signal bold incursions into a new territory. 
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